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ABSTRACT 

An  error  analysis  of  an  hyperbolic  navigation  system  using  high-altitude 
satellites  revealed  that  certain  satellite  arrangements  result  in  long, 
narrow  corridors  on  the  earth  (singular  regions),  in  which  navigation 
errors  are  very  large.  The  study  showed  that  in  these  singular  regions 
the  satellite-user  geometry  results  in  navigation  equations  which  are 
sensitive  to  measurement  errors  and  thus  cannot  be  solved  accurately 
for  all  three  user  coordinates. 
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NAVIGATION  WITH  HIGH-ALTITUDE  SATELLITES: 

A  STUDY  OF  THE  EFFECTS  OF  SATELLITE-USER  GEOMETRY 
ON  POSITION  ACCURACY 


I.  INTRODUCTION 

It  has  become  apparent  that  high-altitude  satellites  have  possibilities  for  navigation  as  well 
as  communication.  Although  use  of  such  satellites  can  allow  the  user  to  obtain  an  accurate  po¬ 
sition  fix,  Schweppe*  pointed  out  that,  for  certain  arrangements  of  the  satellites  and  the  user, 
the  error  in  the  user's  determination  of  his  position  in  three  coordinates  is  extremely  large. 

It  is  the  purpose  of  this  report  to  determine  the  extent  of  such  singular  regions,  i.e.,  regions 
in  which  a  user  cannot  make  a  unique  three -coordinate  fix,  and  to  explore  the  reasons  for  their 
existence. 

II.  SYSTEM  CONCEPT 

The  navigation  system  considered  is  an  hyperbolic  system  in  which  a  user  receives  timing 
signals  transmitted  from  each  of  three  synchronous -altitude  satellites.  The  satellite  clocks, 
from  which  the  timing  signals  are  derived,  are  assumed  to  be  perfectly  synchronized  to  a  master 
clock.  The  user  measures  his  height  as  well  as  the  time  of  arrival  of  each  of  the  satellite  signals. 
In  the  noiseless  case  (that  is,  the  user  has  made  all  his  measurements  perfectly  and  knows  the 
satellite  positions  exactly),  by  using  the  difference  in  the  time  of  arrival  of  the  timing  signals 
emitted  by  two  of  these  satellites,  the  user  can  position  himself  on  the  locus  of  points  satisfying 
this  condition,  that  is,  on  an  hyperbola  of  revolution.  By  taking  one  of  these  satellites  and  a 
third,  he  can  determine  another  hyperboloid.  The  intersection  of  these  two  figures  is  not  suf¬ 
ficient  to  determine  the  user's  position  in  three  coordinates.  The  user  knows,  however,  he  is 
on  the  surface  of  a  sphere  with  radius  equal  to  the  radius  of  the  earth  plus  his  measured  height. 

The  intersection  of  the  sphere  and  the  hyperboloids  will  determine  his  position  in  three  coordinates. 

In  practice,  the  user  is  not  able  to  make  perfect  measurements  and,  therefore,  the  effects 
of  measurement  errors  must  be  considered. 

III.  ERROR  ANALYSIS 

The  quantities  observed  by  the  user  are  (a)  the  time  of  arrival  of  timing  signals  from  each 
of  the  three  satellites  (as  indicated  on  the  user's  clock)  and  (b)  the  user's  height  or  distance  from 
the  center  of  the  earth.  Denote  the  user's  position  in  three-dimensional  space  (origin  at  the 
center  of  the  earth)  by  p  ,  a  three-dimensional  column  vector  (Fig.  1).  Similarly,  denote  the 
satellite  positions  by  the  vectors  p^,  p2,  p^.  If  the  user's  clock  were  synchronized  to  the  master 


*  F.  C.  Schweppe  (private  communication). 
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clock  governing  the  satellite  transmissions  and  there  were  no  errors  in  the  user’s  measurements, 
the  time  of  arrival  of  a  signal  from  a  satellite  would  provide  the  user  with  the  distance  between 
himself  and  the  satellite.  That  is,  the  user  would  observe  indirectly 

VlPu-Ejl  •  j  =  1-2-3  • 

If  the  user  clock  were  offset  from  the  master  clock  by  an  unknown  constant  t  ,  the  user 
would  really  observe 

|p  — P-|+t  =  s.  +  t  ,  j  =  1,  2,  3 

1  £u  -j 1  o  j  o  J 

where  tq  has  the  dimension  of  length,  since  it  has  the  effect  of  changing  the  apparent  distance 
between  the  user  and  the  satellite.  The  user  height  is 

hu=  IpJ 


and  in  the  absence  of  any  noise,  the  observables  are 


y  = 


S.  +  T 
1  O 


S~  +  T 

2  o 


S0  +  T 

3  o 


(1) 


The  quantities  that  y  depends  on  then  are  pu<  p^,  p^,  and  tq.  We  express  these 

quantities  as  a  column  vector,  called  the  state  vector  x, 


£u 


Ei 

x  =  P2 


£3 


(2) 


y  =  h(x)  = 


no  noise. 

h1(x)" 

i  £u  “  £i '  +  To 

h2(x) 

l£u“£zl  +To 

h3(x) 

'  £u  -  £3 1  +To 

h4(x) 

l£u' 

(3) 


The  observables  are,  however,  obscured  by  additive  noise  which  can  be  expressed  in  vec¬ 
tor  form  as 


r  =  y  +  n  =  h(x)  +  n 


(4) 
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We  assume  that  the  a  priori  knowledge  of  the  components  of  the  state  vector  is  in  the  form 
of  a  probability  density.  In  particular,  we  take  x  to  have  statistically  independent  Gaussian 
components  with  mean 

E[x]  =  x 
l-J  -m 


and  covariance 

W  =  E[(x-x  )  (x  —  x  )'] 

where  (  )'  denotes  matrix  transpose  and  E  denotes  the  statistical  expectation  *  The  noise  vec¬ 
tor  n  is  also  taken  to  be  Gaussian  with  zero  mean.  Under  these  assumptions,  the  problem  of 
calculating  the  user  position  from  noisy  observables  can  be  viewed  as  forming  a  statistical  esti¬ 
mate  of  the  state  vector  x.  In  particular,  the  maximum  a  posteriori  probability  estimate  of  the 
state  vector  is  formed,  i.e.,  the  estimate  of  x  that  maximizes  p(x|r).  Let  x  denote  this  esti¬ 
mate.  The  error  analysis  of  such  a  navigation  scheme  consists  of  calculating  the  covariance 
matrix  of  the  error  vector  e  =  (x  —  x). 

In  order  to  do  this,  the  vector  function  h(x)  is  first  linearized  about  the  user’s  true  position 
to  make  this  a  linear  estimation  problem.  Let  the  matrix  H  denote  the  linear  transformation 
required  to  obtain  the  linearized  version  of  the  general  nonlinear  function  h(«),  i.e., 

h(x)  «  h(xQ)  +  H(x  -  xQ)  (5) 

where 


H  =  (h..)  = 


(x=x0) 


This  linear  approximation  of  the  nonlinear  function  h(*)  is  accurate  when  the  difference  between 
the  estimate  x  and  xq  is  small,  as  it  should  be  for  an  accurate  navigation  system. 

The  resulting  error  covariance  matrix 


V  =  (v.j)  =  E[ee'] 
has  been  calculated  as^ 


v  =  W  -  WH[N  +  HWH']"1  HW  =  [H'N-1  H  +  W-1)"1  (6) 

where  W  is  as  previously  defined  and  N  is  the  covariance  matrix  of  the  noise  vector,  i.e., 

N  =  (n. .)  =  E[n.n.) 
ij  1  i  Jj 

where  n^  is  an  element  of  the  noise  vector  n. 

Since  we  are  interested  in  the  estimate  of  p^,  only  a  submatrix  of  V  is  pertinent  to  the 
actual  navigation  errors.  The  upper  left  3x3  submatrix  of  V  is  the  error  covariance  matrix 
of  the  components  of  the  position  estimate  pu-  As  our  measure  of  the  accuracy  of  navigation, 
we  adopt  the  root-mean-squared  (RMS)  error  between  the  estimate  pu  and  the  true  position  pu> 


*  The  statistical  expectation  of  a  matrix  is  taken  to  be  the  matrix  of  expected  values  of  each  of 
the  elements. 

tD.  L.  Snyder,  "Navigation  with  High-Altitude  Satellites:  A  Study  of  the  Errors  in  Position 
Determination,"  Technical  Note  1967-11,  Lincoln  Laboratory,  M.I.T.  (6  February  196*7), 

DDC  648828,  H-803. 
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A  computer  program  was  written  to  calculate  the  V  matrix  and  from  it  the  RMS  position 
error  for  any  positions  of  the  satellites  and  the  user.  The  input  to  the  program  consists  of  the 
a  priori  standard  deviation  of  the  components  of  the  state  vector  x  and  the  observation  noise  n. 
These  are  taken  to  be 


Ranging  error 
Height  error 
User  position  error 
Satellite  position  error 
Clock  error 


10  meters  (RMS) 
10  meters  (RMS) 
10^  meters  (RMS) 
10^  meters  (RMS) 
106  meters  (RMS) 


where  the  square  of  the  ranging  error  and  the  square  of  the  height  error  are  the  diagonal  ele¬ 
ments  of  the  covariance  matrix  N,  the  other  elements  of  this  matrix  being  zero,  since  the  ob¬ 
servation  errors  are  taken  to  be  statistically  independent.  The  squares  of  the  remaining  three 
parameters  are  the  diagonal  elements  of  the  W  matrix,  which  is  also  diagonal  for  the  same 
reason.  In  order  to  eliminate  any  a  priori  knowledge  of  user  position  and  clock  error,  the  var¬ 
iances  of  these  quantities  should  approach  infinity.  This  limit  was  well  approximated  by  using 
10^  meters  as  the  a  priori  variance  of  these  quantities,  which  proved  to  be  of  inestimable  value 
in  view  of  the  finite  word  length  of  the  digital  computer  used  for  the  calculations.  The  other 
parameters  lead  to  reasonable  navigation  accuracies  and  are  realizable  (hopefully)  with  practical 
equipment*  The  program  calculates  the  RMS  error  between  the  maximum  a  posteriori  estimate 
of  pu  and  the  actual  user's  position  for  any  user-satellite  geometries. 

Consider  the  case  where  th*e  satellite  positions  are  as  follows: 


Latitude 

Satellite  1  0.0 

Satellite  2  0.0 

Satellite  3  0.0 


Longitude 

30.0 

0.0 

-30.0 


that  is,  when  all  three  satellites  are  along  the  equator.  The  RMS  errors  for  the  user  at  var¬ 
ious  positions  between  60°  latitude  and  the  equator,  and  40°  longitude  and  0°  longitude  are  given 

2 

in  Fig.  2.  At  60°  latitude,  the  RMS  error  is  of  the  order  of  10  meters  for  all  longitudes  and 
increases  steadily  as  the  user  approaches  the  equator;  at  the  equator,  it  has  a  value  of  10^ 
meters,  but  this  is  the  a  priori  assumption  of  the  position  error  (Fig.  3).  The  observations  and 
calculations  have  not  given  the  user  any  information  about  his  position.  This  is  a  singular  re¬ 
gion  of  the  type  discovered  by  Schweppe. 

We  have  seen  that  the  RMS  error  varies  markedly  over  the  different  user  positions.  The 
question  then  arises  as  to  why  certain  arrangements  of  satellites  and  user  cause  large  errors 
in  position  determination. 


IV.  ANGLE  BETWEEN  HYPERBOLAS 

Consider  again  how  a  user  would  determine  his  position  if  this  were  a  deterministic  case 
(i.e.,  no  noise  and  no  errors  in  satellite  position  knowledge).  He  would  calculate  two  hyperbolas 
of  revolution  and  a  sphere,  all  three  of  which  he  knows  he  must  be  on  (Fig.  4).  The  intersection 
of  these  figures,  therefore,  provides  his  position  fix.  Suppose,  however,  instead  of  being  a 


*  T.  J.  Goblick,  Jr.,  "Navigation  with  High-Altitude  Satellites:  A  Study  of  Ranging  Errors," 
Technical  Note  1966-46,  Lincoln  Laboratory,  M.I.T.  (26  August  1966),  DDC  643851,  H-752. 
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deterministic  case,  there  is  some  error  in  the  user's  measurements  and  thus  an  error  p  in  his 
calculations  of  the  position  of  one  hyperbola.  This  error  would  produce  an  error  in  position 
determination 


€ 


P 

sin  0 


(7) 


where  0  denotes  the  angle  formed  by  the  hyperboloids  intersected  with  the  plane  of  the  earth  at 
the  point  of  the  user.  If  0  is  close  to  90°,  sin 0  is  close  to  one  and  the  error  e  in  the  user's 
determination  of  his  position  is  nearly  equal  to  the  error  p  in  his  calculations.  However,  as 
0  -►  0,  sin  0  -*  0  and  e  -*■  «>,  which  indicates  that  the  error  p  would  be  greatly  magnified  for  small 
angles.  The  hyperboloids  and  the  sphere  intersect  in  such  a  way  that  the  user,  unable  to  meas¬ 
ure  perfectly,  cannot  determine  his  position  accurately  in  three  coordinates.  The  error  in  his 
position  determination  would  be  entirely  dependent  on  his  a  priori  knowledge. 

To  calculate  0,  let  p^  be  a  three-dimensional  vector  (as  in  Sec.  Ill),  origin  at  the  center 
of  the  earth,  denoting  the  user's  position  with  coordinates  (x  ,y  ,  z  ),  and  let  p.,  p?,  and  p- 
represent  the  three  satellites  with  coordinates  (x^y^z^),  z^),  and  (x^y^z^),  respec¬ 

tively.  Then  the  distance  between  the  user  and  each  of  the  satellites  is 

I  Eu  "  Pi  I  =  n/(Xu  "  Xi)2  +  (yu  "  ^i)2  +  <zu  “  Zi)2  =  Si  .  i=  1.2,3  .  (8) 

Since  the  time  required  for  the  user  to  receive  the  timing  signals  is  a  function  of  the  dis¬ 
tance  between  the  user  and  the  satellite,  we  can  describe  the  hyperboloid  on  which  the  user  is 
located  as  the  locus  of  points  such  that  the  difference  in  the  distance  between  that  point  and  two 
other  points  is  constant.  One  hyperboloid  is  determined  by  satellite  1,  satellite  2,  and  |  s^  —  s^l; 
and  another  by  satellite  2,  satellite  3,  and  |  s^  —  s^) .  From  Eq.  (8),  the  two  hyperboloids  are 

f l  (X,  y,  z)  =  I  sl  -  S2 1  =  J(x  -  x^2  +  (y  -  y1)2  +  (z  -  z^2 

-  J(x  -  x2)2  +  (y  -  y2)2  +  (z  -  z2)2  (9) 

f2(x,  y,  z)  =  |  s2  -  s3 1  =  J(x  -  x2)2  +  (y  -  y2)2  +  (z  -  z2)2 

-  J(x  -  x3)2  +  (y  -  y3)2  +  (z  -  z3)2  .  (10) 

If  f  is  a  function  of  three  independent  variables  (x,  y,  z),  the  gradient  of  the  function  is 
defined  as 


where  i,  j,  and  k  are  unit  vectors  along  the  positive  x,  y,  and  z  axes,  respectively.  Geomet¬ 
rically,  Vf,  evaluated  at  the  point  (xq,  yQ,  zq),  is  a  vector  whose  direction  is  normal  to  the  level 
surface  [i.e.,  the  set  of  all  points  such  that  f(x,  y,  z)  =  c]  at  the  point  (xq,  yQ,  zq).  In  particular, 
f^(x,y,  z)  in  Eq.  (9)  defines  a  family  of  hyperboloids,  the  value  of  |  s ^  —  s2 1  defining  the  particu¬ 
lar  hyperboloid.  The  gradient  of  f^  evaluated  at  a  point  (xq,  y  ,  zq)  will  be  a  vector  whose  direc¬ 
tion  is  normal  to  the  particular  hyperboloid  passing  through  (xq,  yQ,  zq).  Likewise,  we  could 
find  a  vector  whose  direction  is  normal  to  the  hyperboloid  defined  by  f2(x,  y,  z)  in  Eq.  (10)  pass¬ 
ing  through  the  point  (x  ,  yQ,  zq). 

If  Vf  and  Vf2  are  the  gradients  of  the  functions  f^  and  f^,  then  the  angle  0,  formed  by  the 
intersections  of  the  hyperboloids  and  the  tangent  plane  of  the  sphere,  can  be  found  by  computing 
the  angle  projected  in  the  tangent  plane  between  Vf^  and  Vf^  evaluated  at  the  user’s  position. 


Let  E  be  a  unit  vector  normal  to  the  plane  tangent  to  the  sphere  at  the  point  of  the  user.  Then 
Vq  X  E 

M1  =  I  Vf*  xe| 

and 

Vf2  x  E 
-2  =  |  Vf2  x  E| 

are  unit  vectors  in  the  plane  tangent  to  the  sphere,  and 
cos  0  =  u^  •  u^ 

or 

0  =  cos"1  (u  •  u2)  (13) 

When  Eq.  (13)  is  evaluated  for  satellite -user  geometries  where  navigation  accuracy  is  known 
to  be  poor,  the  values  of  0  should  be  near  or  exactly  zero.  In  particular,  if  this  expression  is 
evaluated  for  the  satellite  geometry  used  in  Sec. Ill  at  the  user  positions  that  had  the  maximum 
RMS  position  error,  that  is,  when  the  user  is  located  at  any  point  along  the  equator,  cos  0  does 
equal  one  or  the  angle  between  the  hyperbolas  in  the  plane  tangent  to  the  earth  equals  zero  (see 
Fig.  5). 

Since  it  is  more  difficult  to  solve  this  equation  for  non -equatorial  cases,  a  computer 
program  was  written.  Input  to  the  program  consists  of  the  longitude  and  latitude  of  the  three 


(11) 


(12) 


6 


satellites  and  the  user;  l/sin9  is  the  output  quantity  rather  than  9,  since  the  error  magnifica¬ 
tion  is  seen  in  Eq.  (7)  to  depend  on  l/sin9. 

There  is  a  direct  correspondence  between  tl\e  user  locations  resulting  in  large  values  of 
l/sin9  and  user  locations  with  large  RMS  position  errors  for  the  same  satellite  geometry.  The 
expression  l/sin9  had  a  minimum  value  at  60°  latitude  for  all  longitudes  and  increased  steadily 
until,  along  the  equator,  9  became  zero  and  l/sin9  could  not  be  computed.  It  was  also  at  60° 
latitude  that  the  minimum  RMS  position  error  occurred,  and  along  the  equator  that  the  maximum 
occurred. 


V.  RESULTS 


It  is  interesting  to  explore  satellite-user  geometries  with  less  symmetry  than  the  equatorial 
case  of  Fig.  5.  The  two  computer  programs  were  run,  therefore,  for  other  satellite  and  user 
positions. 

The  following  satellite  positions  provide  a  case  with  only  east-west  symmetry. 


Satellite  1 
Satellite  2 
Satellite  3 


Latitude 

Longitude 

0.0 

30.0 

10.0 

0.0 

0.0 

-30.0 

Navigation  accuracy  for  user  positions  from  +60°  to  —60°  latitude  and  from  —40°  to  +40°  longi¬ 
tude  was  studied.  A  singular  region  again  appeared.  The  maximum  RMS  error  of  the  user  at 
—  40°  longitude  occurred  at  latitude  +41°;  at  0°  longitude,  it  appeared  at  49°  latitude;  and  at  40° 
longitude,  the  maximum  was  at  41°  latitude.  The  effect  of  moving  the  satellite  north  was  to 
move  the  singular  region  in  that  direction  while  the  southern  hemisphere  became  free  from  any 
singular  points.  The  actual  values  of  the  errors  are  meaningful  only  in  a  relative  sense  because 
the  region  was  not  sampled  finely  enough  to  conclude  that  the  worst  point  was  found.  The  results 
of  the  program  to  calculate  the  angle  between  the  hyperboloids  for  this  same  geometry  again 
showed  the  maximum  value  of  l/sin9  occurring  at  exactly  the  same  user  positions  as  the  maxi¬ 
mum  RMS  position  errors  [see  Fig.  6(a-b)]. 

If  satellite  2  is  now  moved  to  15°  latitude,  15°  longitude,  all  symmetry  is  destroyed,  and 
user  positions  from  -40°  to  +40°  longitude  and  +75°  to  —75°  latitude  must  be  studied.  (The  in¬ 
crease  in  the  range  of  latitude  is  needed  because  of  the  more  northerly  position  of  the  satellite.) 
Again,  there  is  a  narrow  locus  of  singular  points  with  the  worst  point  this  time  occurring  (user 
longitude  0°)  at  user  latitude  68°.  This,  however,  is  nearly  out  of  the  region  of  mutual  visibility* 
of  all  three  satellites  [see  Fig.  7(a-b)].  The  region  of  singular  points  could  probably  be  moved 
entirely  out  of  the  region  where  the  user  could  see  all  three  satellites  by  further  adjustment  of 
the  northern  satellite. 

We  can  conclude  that  the  region  of  singular  points  exists  whether  the  satellites  are  arranged 
symmetrically  or  unsymmetrically.  The  arrangement  of  the  satellites  does,  however,  have  an 
effect  on  the  location  of  the  region  of  large  errors. 


*  In  this  report,  visibility  curves  are  computed  using  an  elevation  angle  of  7°. 
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Taking  the  positions  of  the  satellites  to  be 


Satellite  1 
Satellite  2 
Satellite  3 


Latitude 

0.0 

0.0 

0.0 


Longitude 

45.0 

0.0 

-45.0 


we  find  that  the  singular  points  again  occur  on  the  equator  as  expected,  and  that  the  maximum 
values  of  the  RMS  position  errors  are  the  same  along  the  equator  as  for  the  previous  equatorial 
case  because  of  the  a  priori  assumptions.  If,  however,  we  compare  the  RMS  values  of  the  user 
positions  off  the  equator  for  the  30°  and  45°  cases,  we  find  that  the  corresponding  values  for  the 
3  0°  case  are  2.1  times  larger  than  those  for  the  45°  case,  and  that  l/sin0  is  0.5  times  larger 
for  the  30°  case  [see  Fig.  8(a-b)].  The  effect  of  having  the  equatorial  satellites  90°  instead  of  60° 
apart  was  a  decrease  in  the  width  of  the  region  of  large  errors.  However,  the  region  of  mutual 
visibility  when  the  two  satellites  are  90°  apart  is  so  limited  that  three  satellites  would  not  be 
sufficient  to  cover  the  North  Atlantic,  for  example;  whereas,  this  could  be  done  if  the  satellites 
were  only  60°  apart. 

Thus,  if  high-altitude  satellites  are  to  be  used  for  navigation,  consideration  must  be  given 
to  arranging  the  satellites  so  that  the  user  will  not  be  in  a  singular  region,  that  is,  so  that  the 
user  will  be  able  to  determine  his  position  accurately.  But  consideration  must  also  be  given  to 
making  the  area  of  mutual  visibility  of  the  user  and  the  satellites  large  enough  for  the  system  to 
be  practical. 


VI.  SUMMARY 

From  the  preceding,  the  following  conclusions  can  be  drawn: 

(1)  It  appears  that  all  satellite  geometries  result  in  long,  narrow  corridors 
in  which  navigation  accuracy  is  poor. 

(2)  The  cause  of  this  singular  region  is  geometric  in  nature;  that  is,  the  two 
hyperboloids  and  the  sphere  which  locate  the  user  intersect  in  such  a  way 
as  to  greatly  magnify  any  error  in  the  user’s  measurements. 

(3)  A  system  with  all  satellites  on  the  equator  always  yields  a  singular  region 
along  the  equator.  Thus  a  practical  system  would  utilize  at  least  one  non- 
equatorial  satellite. 

(4)  If  two  of  the  satellites  are  on  the  equator  and  the  third  is  located  far  enough 
north,  the  singular  region  can  be  shifted  outside  the  range  of  mutual  visi¬ 
bility. 

(5)  If  two  of  the  satellites  are  on  the  equator,  the  singular  region  is  located  in 
the  same  hemisphere  as  the  third  satellite;  i.e.,  the  opposite  hemisphere 
will  be  free  from  singular  points. 
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Fig.  2.  RMS  position  errors  for  all  equatorial  satellites.  (All  values  in  meters.) 


Fig.  3.  Cross  section  of  RMS  position  error  vs  latitude  (longitude  =  0°). 
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Fig.  4.  Intersection  of  hyperboloid  with  tangent  plane  of  earth. 
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Fig.  5.  1/sin 0  for  all  equatorial  satellites 
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Fig.  6(b).  RMS  position  errors  for  an  east-west  symmetric  satellite  geometry. 
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